Dependence of the superconducting critical temperature on the number of layers 

in homologous series of high-T^ cuprates 



O 

o 

(N 

in 

o 
o 

u 



T. A. Zaleski and T. K. Kopec 

Institute of Low Temperature and Structure Research, 

Polish Academy of Sciences 

P.O. Box 14.10, 50-950 Wroclaw, Poland 

We study a model of n-layer high-temperature cuprates of homologous series like 
HgBa2Can-iCun02+2n+<5 to explain the dependence of the critical temperature Tc (n) on the number 
n of Cu-O planes in the elementary cell. Focusing on the description of the high-temperature super- 
conducting system in terms of the collective phase variables, we have considered a semi-microscopic 
anisotropic three-dimensional vector XY model of stacked copper-oxide layers with adjustable pa- 
rameters representing microscopic in-plane and out-of-plane phase stiffnesses. The model captures 
the layered composition along c-axis of homologous series and goes beyond the phenomenological 
Lawrence-Doniach model for layered superconductors. Implementing the spherical closure relation 
for vector variables we have solved the phase XY model exactly with the help of transfer matrix 
method and calculated Tc (n) for arbitrary block size n, elucidating the role of the c-axis anisotropy 
and its influence on the critical temperature. Furthermore, we accommodate inhomogeneous charge 
distribution among planes characterized by the charge imbalance coefficient R being the function of 
number of layers n. By making a physically justified assumption regarding the doping dependence 
of the microscopic phase stiffnesses, we have calculated the values of parameter i? as a function of 
block size n in good agreement with the nuclear magnetic resonance data of carrier distribution in 
multilayered high-Tc cuprates. 
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I. INTRODUCTION 

Characteristic common structure of high-temperature 
superconducting cuprates has been puzzling researchers 
since the discovery of the phenomenon by Bednorz and 
Miilleri It is now accepted that all relevant carriers of 
both spin and electricity derive from the hybridized an- 
tibonding orbitals in the Cu-O planesi^ Obviously, the 
interlayer structure cannot be ignored^ In the normal 
state, due to huge anisotropy, coherent transport in 
the direction perpendicular to Cu-O planes is blocked. 
However, it is observed that Tc is strongly dependent 
on the interlayer structure, which means that three- 
dimensional (3D) coupling of planes must play an im- 
portant role in the onset of superconductivityj^ It sug- 
gested to apply the Lawrence and Doniach model (LD) 
to cuprates, which originally was proposed for layered su- 
perconductors and is based on the Ginzburg-Landau de- 
scription with neighboring layers coupled by the Joseph- 
son tunneling^ This approach was pursued to study 
Josephson plasma resonance?^ anisotropy of critical fields 
Hci and Hc2r^ vortex matter behavior iS However, the 
biggest disadvantage of the LD model is that it is based 
on mean-field approximation, which neglects the role of 
phase fluctuations, the latter playing a profound role in 
cupratesiiS 

There is a significant evidence that in cuprates the 
Cu-O planes are Josephson coupled and the anisotropy 
of c-axis to a6-plancs penetration depth can vary from 
5 to even 250 in various high-Tc compounds. This gave 
rise to interplane theories e.g. interlayer tunneling model 
(ILT) by Andersoufii which assumed that carriers were 
incoherent between layers in the normal phase, how- 
ever coherent transport of Cooper pairs in a form of 



Josephson tunneling in superconducting state was al- 
lowed. Unfortunately, the original version of the ILT 
theory by Anderson is no longer considered as a viable 
model of high-Tc materials, since it provides no more than 
1% of the condensation energy in certain cupratesiiSii^ 
However, the Josephson coupling between Cu-O planes 
is still believed to play an important role in cuprate 
oxidesii^ In the process of synthesis of high-Tc mate- 
rials it is possible to control (to some extent) the dis- 
tances between copper-oxide planes or the content of 
the regions that separate them. It allows to study, how 
such changes influence observed properties of these com- 
pounds like critical temperature or electrical resistivity 
and elucidate the role of the inter-layer coupling pi^*i& 
In mono- or bi-layered high-Tc materials (like LSCO or 
YBCO compounds), all copper-oxide planes are equiva- 
lent, which means that charge carriers are distributed 
homogeneously among copper-oxide planes. However, 
multi-layered high-Tc materials are characterized by the 
existence of crystalographically inequivalent outer (OP) 
and inner (IP) Cu-O planes, which give rise to inhomo- 
geneous carrier distribution among them. This, in turn, 
may have a profound influence on the critical tempera- 
ture, which is know to be strongly doping dependent. For 
example, high-temperature superconducting homologous 
series, like HgBa2Ca„_iCu„02+2n-i-<5, are materials shar- 
ing the same charge reservoir, but differing by the num- 
ber n of Cu-O planes in the unit celliiLi^ They also have 
the same characteristic dependence of the critical tem- 
perature as a function of n with a maximum for n — 3 
or n = 4. Although, these groups of cuprates have been 
extensively studied, the origin of this extremum is still 
not fully understoodii2i22iSi 

Recently, Chakravarty et al. have argued that the ex- 



istcnce of the competing order parameter that is nucle- 
ated by charge imbalance is ultimately responsible for the 
critical temperature downturn with m^ They present a 
mean-field theory for a single block of n superconduct- 
ing layers, however neglecting the inter-block coupling. 
Moreover, they use phenomenological Ginzburg-Landau 
theory at zero-temperature to anticipate the values of the 
critical temperature Tc (n) , based on the magnitude of 
the T = superconducting order parameter. Charge im- 
balance is introduced as a competing order parameter, 
e.g. d-density wave^ It depresses the critical temper- 
ature for n > 2, thus experimentally observed shape of 
Tc (n) is obtained. However, as the authors note, the 
phase fluctuations, which are not taken into account, 
might depress Tc for n = 1 even more than for n = 3. 
Although, it should correct the results for 1 < n < 3, it 
is not obvious, whether the maximum of Tc (n) can be 
observed at all. 

It is well known that binding of electrons into pairs is 
essential in forming the superconducting state, however, 
its remarkable properties - zero resistance and Meissner 
effect - require the phase coherence among the pairs as 
well. In conventional BCS superconductors the phase or- 
der is unimportant for determining of the value of the 
transition temperature Tc. However, in superconduc- 
tors with low carrier density such as high- Tc oxide su- 
perconductors, phase fluctuations may have a profound 
influence on low temperature properties »iii In particular, 
for cuprate superconductors, the conventional ordering 
of binding and phase stiffness energies appears to be re- 
versed. Thus, the issue how pairing and phase correla- 
tions develop is a central problem of high- Tc supercon- 
ductivity. The measurements of the frequency dependent 
conductivity, in the frequency range 100-600GHz, shows 
that phase correlations indeed persist above Tc, where 
the phase dynamics is governed by the bare microscopic 
phase stiffnesses iSi 

In the present paper we propose a semi-microscopic 
model of n-layer cuprates, which is founded on micro- 
scopic phase stiffnesses that set the characteristic en- 
ergy scales: in-plane J|| , inter-plane (in-block) Jx and 
inter-block Jj_ couplings. We employ an approach that 
goes beyond the mean field level and is able to cap- 
ture both the effects of phase fluctuations and huge c- 
axis anisotropy on the superconducting phase transition. 
We show that the depression of the critical temperature 
for n > 3 can be obtained naturally without invoking 
any competing (hidden) order. Since superconducting 
carriers are located in Cu-0 planes, the in-plane phase 
stiffness must naturally be doping S dependent. As this 
dependence can be simply deduced from the established 
Tc — S phase diagram for cuprates, introduction of charge 
imbalance via the doping dependence of microscopic in- 
plane phase stiffness allows us for explanation of Tc (n) 
as a function of block layers n in perfect agreement with 
experimental findings^ 

The outline of the reminder of the paper is as follows. 
In Section II we analyze a single-layer system. We con- 



struct an anisotropic three-dimensional XY model, which 
will be solved in the spherical approximation to explicate 
the dependence of critical temperature on the anisotropy 
rj = J±/J\\ along c-axis. Subsequently, in Section III we 
elaborate on a model of stacked n-block layers relevant 
for high-temperature superconducting homologous series. 
We map the system onto a spherical XY model, which is 
solved exactly for arbitrary n with the help of the transfer 
matrix method. In Section IV we determine the critical 
temperature Tc (n) as a function of number of layers n. 
In Section V we incorporate the inhomogeneous charge 
distribution in homologous series by introducing doping 
dependent in-plane phase stiffnesses Jy (6) for the sepa- 
rate n layers in the block structure. We found that the 
maximum of the phase coherent superconducting transi- 
tion will be a bell-shaped curve as a function of n with 
a maximum for n — 3 for appropriate values of outer 
{6op) and inner (Sip) planes dopings within a block. Fi- 
nally, in Section VI we summarize the conclusions to be 
drawn from our work. Some supplementary material is 
relegated to the Appendix. 



II. SINGLE LAYER CUPRATES 

In underdoped high-temperature superconductors, two 
temperature scales of short-length pairing correlations 
and long-range superconducting order seem to be well 
separatediS. We consider the situation, in which local 
superconducting pair correlations are established and the 
relevant degrees of freedom are represented by phase fac- 
tors < (fie (vi) < 2tt, where r^ numbers lattice sites 
within i-th ab plane. The system becomes superconduct- 
ing once U(l) symmetry group governing the ipi (ji) fac- 
tors is spontaneously broken and the non-zero value of 
^giVf(ri)^ appears signahng the long-range phase order. 
The dependence of critical temperature on the spatial c- 
axis anisotropy is of paramount importance for cuprate 
systems, in which two-dimensional physics predominate. 
First, we consider a single-layer system described by the 
following Hamiltonian: 

H[ip] = -J||^^cos[(p£(r,)-^£(rj)] 

-Ji-^^cos[ipi{v{)-ipi+i{v,)], (I) 

t i 

where J\\ > and Jj_ > are in-plane and inter-plane 
microscopic phase stiffnesses, respectively. The factors 
fi (ri) are placed in sites of a three-dimensional cubic 
lattice, where i = l,...,iV|| with A'^n being the number 
of sites in a plane, £ ~ I, ..., N±, where N± denotes the 
number of layers and N = A^iiA^^ is the total number of 
sites. The partition function of the system reads: 



Z = 



2ir 



■0Hlv] 



(2) 
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Figure 1: (Color online) Pictorial representation of the var- 
ious crossovers due to increased spatial anisotropy between 
in-plane and inter-plane couplings: 3D^2D (isotropic inter- 
actions to uncoupled planes) and 3D— >1D (isotropic interac- 
tions to uncoupled chains). 



Si (Yi) variables runs only over the values, which satisfy 
the unit length condition S^ (r^) = 1. Unfortunately, the 
partition function in Eq. ij^l cannot be evaluated exactly. 
However, replacing of the rigid length constraint in Eq. 
(0) by a weaker spherical closure relatio n^^i^^ 



5[Sl{v,)-l] 



^Es|(rO-l 



(6) 



renders the model in Eq. Q exactly solvable. The rela- 
tion in Eq. © means that the unit length of the S^ (r.^) 
vectors is maintained on average. The spherical con- 
straint can be conveniently implemented with the help 
of Dirac-(5 function and the partition function in Eq. (0) 
becomes: 



Z = 



+ 00 



X{cfSdv,)5 



i.i 



J2sUr^)-N 



i,£ 



,^PH[S] 



(7) 

Furthermore, the i5-function can be conveniently repre- 
sented in a spectral form by the appropriate integral 
S (x) = /_ d(/2TTicxp{—(x), which introduces the La- 
grange multiplier (^. Consequently, Eq. |(7J) can be writ- 
ten as: 



/■+°° dC 
Z = j ^exp[-7V/3/(C)], 



(8) 



where the free energy per site is given by: 



/(C) 



c 

/3 N(3 



-In/ n^'Sar.) 



X exp<^-C5]S|(r,)-/3i?[S] 



(9) 



where j3 = l/ksT with T being the tempera- 
ture. Introducing two-dimensional vectors S^ (r,) = 
[Sxe (r,;) , Syi (r,;)] of the unit length S^ (r,;) = S^^ (r,) + 
Sy^ivi) = 1 defined by 

Si (r,;) = [cos ^i (ri) , sin ipi (r^)] , (3) 

the Hamiltonian can be written in the XY-model form: 

H[S] = -J||^5]S,(r,;)-S,(r,) 

-•^^EEsarz)-S,+i(r,). (4) 

I i 

In terms of the vector variables the partition function in 
Eq. becomes: 



Z = 



+ 00 



X{{d'Si{v.,)5[SUv.)-l]]e~P"^^\{h) 



i,i 



where d^S^ (r^) = dSj;i{ri) dSyi{ri) and the Dirac-5 
function S [Sj (r^) — l] assures that the integration over 



In the thermodynamic limit N —>■ oo the dominant con- 
tribution to the integral in Eq. ((HJ comes from the saddle 
point C = Co of / (C) with value of Co determined from 
the condition: 



dfiC) 



dC 



0. 



(10) 



C=Co 



To evaluate the free energy of the system f = f (Co) we 
diagonalize the Hamiltonian from Eq. Q) by introducing 
Fourier transform of Si (ri) variables: 



Si (r,) = 



1 



NnNt 



k.q. 



Sk,9,e 



— i{]<.ri-\-qzCi 



(11) 



where k = (k^, ky) with —-< k^, ky < - are wave vec- 
tors within a6-planes, while — - < gz < - is wave vector 
along c-axis. Here, a and c are lattice constant within ab- 
plane and c-axis, respectively. The Hamiltonian in Eq. 
Q then transforms to: 



H = -:5^E'^(k'9^)Sk*-S- 



-9^' 



(12) 



k,gz 



where for the Fourier-transformed interactions Jy and J± 
on the 3D anisotropic cubic lattice with nearest-neighbor 
coupHng we readily find: 

J (k, Qz) = 2 J|| [cos {akx) + cos {aky)] + 2J± cos {cq^) . 

(13) 
As a result, the integral in Eq. Q becomes Gaussian and 
the free energy can be calculated: 

+ CJO 



f3 f3N J_^ li "''= 



X exp ■ 



-y 



k,g. 



1 



d^S 



k,9. 



Co - -Jik,q;, 
Kg. 
dReS^ (iImSk,5^ 



Sk,g^Sk,g^ 



(14) 



The emergence of the critical point in the spherical model 
is signaled by the conditionS^ 



G-i(k = 0,9. = 0) = Co-;^/3Jo = 0, 



(15) 



which fixes the saddle-point value of the Lagrange mul- 
tiplier C within the ordered phase: 



Co - 1/3 Jo, 



(16) 



where Jq — J (k = 0, q^ ~ 0) . Here, G (k, q^) is the order 
parameter susceptibility defined by: 



G (k, q^) = (Sk,g, • S_k,-gJ 
with the statistical average: 

1 f + OO 



(17) 



Z 



xS 



ySJ{r,)~N 



i,£ 



..e 



-PHIS] 



(18) 



From Eqs. (|10|l and p4|l we obtain the critical tempera- 
ture in the form: 



/3c 



E 



1 



lim , 

AT^oo N f^ Jo- J (k, q^ 



(19) 



For one (ID) and two-dimensional (2D) systems with 
short-range interactions the ordered phase exists only in 
the zero temperature limit, hence for dimensions d < 
2 the spherical model docs not exhibit phase transi- 
tion with Tc > 0, in agreement with Mermin-Wagner 
theoremiS^ For d = 3 the critical temperature is finite 



2kBTc 
Jo 



0.659, 



(20) 



where the factor 2 in front of ksTc comes from the num- 
ber of components of the vector XY modcli^ 
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Figure 2: Superconducting critical temperature Tc as a func- 
tion of anisotropy r\ for crossovers from isotropic 3D to ID or 
2D anisotropic systems in logarithmic scale for r\. 



Introducing spatial anisotropy parameter < ?/ < 
1 (being the ratio between J|| and Jj^). J(k, g^) 
can be written in a form to accommodate 3D— >2D 
crossover froirr the isotropic superconductor on the three- 
dimensional cubic lattice (77 = 1) to the uncoupled 
(r\ = 0) copper-oxide planes (see. Fig. |2J|: 

J(k, (72) = 2 J [cos (afca;) + cos (afcy) -f 77 cos (c^z)] 

J=Jj|, v^J^/J\\- (21) 

In a similar way we can study a BD^ID crossover 
from the isotropic superconductor to the uncoupled one- 
dimensional chains: 

J{k,qz) = 2J{r/ [cos(afc^) -f cos(afcy)] -h cos(cgz)} 

J=Jx, '? = J||/J±. (22) 

The critical temperature can be found for both cases from 
Eq. 1191 by performing momentum integration. In the 
thermodynamic limit we have: 



lim — 



E' 

Kg. 



d^k 



-^ i27r/ay J-^ (2^/c) 



dqz 



(23) 



For 3D^2D crossover it reads: 

Jo _ 2 (2 + T]) {2y/TTv - V^T2^) 



ksTc (77) 
xK 



777r^ 



2^T+7j - ^I+2rj) v/4 + 2^+2 



277 

^ (2 VTT^ - V'4 + 2^) ( a/4 - 



2r? 



(24) 



while for 3D^1D changeover the Eq. (|19() : 



Jo 



ksTc (rj) 



2 (1 + 2i]) {2^TTv - V2~4^) 
7r2^ 



xK 



xK 



2lV^-\/2 



yi^i-vTT^j 



2{Vv 




v^ 



(25) 



where K (x) is the complete elliptic integral of the first 
kind:-- It is interesting, how these Tc {rf) values behave 
in the limiting cases. For the 3D^2D crossover from Eq. 
(1^ we obtain: 



Jo 



1.319 - 0.094477^ - 0.0559?/-^ for?/ -^ 1 

27r TT-q [In (82/77) + 2] 



In (32/,;) 



2 In (32/?7) 



for 77 ^ 0, 



(26) 
while for 3D^1D changeover, Eq. (|25|) it follows that: 



ksTc 



1.319-0.0944772 

{2 + V2)tt^^ 

4K2 (3 - 2\/2) 



O.IO6774 for?7^ 1 
for ?7 -^ 0. 



(27) 



The results are presented in Fig. |21 and in Table Q] 
When the in-plane interactions are reduced, the criti- 
cal temperature drops quite rapidly as a function of 77 
compared to the isotropic system and for 77 = 10~^ is 
less than 1% of the isotropic value. However, in the 
case of the reduced coupling along c-axis (weakly cou- 
pled planes), even for anisotropy parameter as small as 
Tj = 10^® the critical temperature is still about 30% of 
the isotropic three-dimensional value. This result is very 
important as it proves that even very small inter-planar 
coupling (sometimes too small to be directly observed ex- 
perimentally) renders the phase transition in the 3D uni- 
versality class ^ with the critical temperature decreas- 
ing slowly like Tc ~ 1/ log 77. This result raises serious 
doubts whether the outcome of intercalation experiments 
in cupratesji^ which show small changes of the critical 
temperature with increase of the inter-layer distance can 
be interpreted as an argument in favor for a purely 2D 
nature of high- Tc superconductivity. 



III. MULTI LAYERED SYSTEM 

To describe superconducting cuprates of homologous 
famihes, we modify the model presented in Sec. E] to 
include the presence of the n-layer block structure. To 
this end, we group every n consecutive layers in blocks. 
The intra-layer (in the ab planes) interactions are given 
by J|| > 0, while the plane coupling along c direction 
is different inside J± > and between J , > blocks. 
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Table I: Ratio of the superconducting critical temperature 
Tc (77) to the isotropic three-dimensional value Tc (1) for sev- 
eral values of the anisotropy parameter 77 facilitating the 
crossovers from isotropic 3D to anisotropic ID or 2D systems. 



respectively (see, Fig. (S)). The phase Hamiltonian of the 
system can now be written: 



H = - Jjj Y^ Y^ cos [ifi (r,) - ifi (rj)] 
i i<j 

n-2 

-J± X! X! X! *^°^ [frnn+a {^1) - ^„in+a+l (r^)] 
i m Q>0 

- J± ^ Yj '^'^^ [Vmn+n-1 (fj) - 'fmn+n (r^)] , (28) 
i m 

where indices i and j run over sites in the a6-plancs of 
the £-th layer, while a numbers the layers within the 777- 
th block. Introducing two-dimensional vectors from Eq. 
Q, the Hamiltonian reads: 

^ = -j|iEEs^(^'')s^(^j-) 

n-2 

J± / J / ^ / ^ ^mn+a {^i) ^mn+a + 1 [^i) 
i m Q>0 

— J_L 2^ A^ ^inn+n~l (l";) S„i„+„ (ri) . (29) 



Furthermore, we make use of the spherical closure re- 
lation in Eq. ©. However, the direct implementation of 
the Eq. ^ is obstructed by the lack of complete transla- 
tional symmetry: grouping of layers within blocks breaks 
translational symmetry along c axis, since the inter-plane 
coupling varies with period of 77, when moving from one 
plane to another. As a result, the three-dimensional 
Fourier transform of the S^ (r^) variables in Eq. Hll|l . 
used to diagonalize single-layer Hamiltonian in the previ- 
ous Section, cannot be performed. To overcome this dif- 
ficulty, we implement a combination of two-dimensional 
Fourier transform for in-plane vector variables 



1 



II k 



(30) 



and transfer matrix method for one-dimensional deco- 
rated structure along c-axis^ This operation diagonal- 
izcs all terms in the Hamiltonian in Eq. H29|l with respect 




to k, leaving the dependence on layer index (, unchanged. 
Therefore, the partition function can be written in the 
form: 



Z = 



-foo 



^exp<JA.C+^ln 



+ 00 



n^'Sk^ 



X exp 



-]^ESk.^lWSk.' 



(31) 



Figure 3: (Color online) Structure of the multi-layer 
high-temperature superconductors of homologous series of 
cuprates. Coupled Cu-O planes are grouped in blocks of n 
planes. The microscopic phase stiffnesses Jy , J_l and J^ are 
indicated by the thick lines. 



where A|^ (k) is an element of a square N±^ x N±_ band 
matrix, appearing as a result of non-trivial coupling 
structure along c-direction: 



J 



An^ (k) 
































>^i 



(32) 



N^ 



where 

J|| (k) = 2 J|| [cos (afc^) + cos {aky)] . 
Here, Anj_ has square n x n diagonal blocks: 



(33) 



B„ 



c 



/3J||(k) 
2 

2 






Ma 

2 







_I3J± 



C- 






;3J||(k) 



containing interactions inside the stack of n planes. Per- 
forming Gaussian integration in Eq. 131|l . one arrives at 
the equation involving determinant of the matrix in Eq. 



Z = 



2TTi 



exp 



iVC + ^ln 



N, 



,N, 



del A 



N, 



(35) 



Thus, the problem of finding the partition function in Eq. 
H35(l is to find a way to calculate the determinant of the 
matrix An '■ 



det(^w^) = |^ 



Ni 



(36) 



However, because of the band structure of the ma- 
trix An^ it can be conveniently done using generalised 
Laplace method. As a result, the determinant in Eq. H3ti|) 
can be written in a recurrence form: 



1^ 



N, 



\Bn\ \An^^ 



(3Ji_ 



l-Sn-ll |^Afj_-r 



(37) 
where A^^^n and A^^-n-i arise from removing the first 
n and n + 1 rows and columns of A^^ , respectively. Sim- 
ilarly to Eq. H37|) . for |A7v^_i| the recurrence formula 
can be written as: 

I^AT^-ll = l-Sn-ll |^Afj_-n| 

[2) \Bn-2\\AN^-n-l\. (38) 

Defining |i?o| = 1, the formula in Eq. (|38|l holds for any 
n > 2. The results for the single-layer system (77, = 1) 
from the previous Section can be recovered by performing 
calculations for any n > 2 and setting J± = J^. Using 
Eqs. l|?7|l and l|5S|l we can write the recurrence relation 
for the Ajs[^ in Eq. H32() as an operator equation: 



\AnA 
\An^-i\ 



|^Afj_-n| 
I^AT, -n-l| 



(39) 



where Ti is transfer matrix: 



Ti 



\Bn\ 
\Bn-l\ 



PJ± 



[iJ. 



\Bn- 



\Bn-2 



(34) Applying the rule in Eq. 



(40) 



recursively, one obtains: 






--N, -1 



|A„_i| 



(41) 



Therefore, finding the solution of the Eq. (|41|l reduces 
to the problem of determining the power of the transfer 
matrix operator Ti. To this end we diagonalize the ma- 
trix Ti with the help of a unitary transformation Ui (see. 
Appendix ^ : 



Ti=Ui 



Ai 
A2 



^r 



(42) 



where Ai and A2 are eigenvalues of Ti, whereas Ui is a 
unitary matrix built of eigenvectors of Ti. Therefore, 



Tf''--'=U, 



A{ 



^N, -1 



^N^ 



U^\ 



(43) 



The eigenvalue Ai is chosen to satisfy Ai > A2. In the 
large-A'^x (thermodynamic) limit only the largest eigen- 
value Ai will contribute to the statistical sum in Eq. H35|l . 



Ai 
A2 



J-Af^-l 



Ai 



= K 



N± 



1 
A2/A1 

1 




iAT^-l 



(44) 



Combining Eqs. |(1T|) . l(^ and (|^ we obtain finally: 



\AnA 



N±~*OD 



A?"--V 



1 




^r 






(45) 



where: 




/3J, 



\Bn-2\ + 



\ 



\B„ 



pj'l_ 



\Br,n 



(/3Jl)'|S„_i| 



(46) 



From Eq. (|45|) one can extract the value of determinant of the matrix Apf^ : 

--N, -1 



\AnA = ^I ^ (^)n [(^r')n \Bn\ + {U^'),, |i?„-i|] 



(47) 



where {U) ^i, and {lAi ) denote elements of the matrix Ui and the inverse matrix Ui , respectively (see, Appendix 
Ej). Explicitly, 



\A 



N, 



\BJ - 



PJi 



\B, 



ri-2 




i7Vi-l 



(48) 



is the function of determinants of matrices i3„ with intra- 
block interactions J^ and 

C = {U,)^, [(Wr')n \Bn\ + K')2i \Bn^i\] ■ (49) 

The determinant of the matrix i3„ in Eq. H34|l can be 
found using similar method by employing the Laplace ex- 
pansion (see Eq. H37|) - (|45|l 'l. However, since the order of 
Bn is finite, both eigenvalues of the new transfer matrix 
T2 must be taken into account: 



l-S„+i| 

\Bn\ 



1/2-1 



1/2-1 
f 

A^ 



\B2\ 
\Bi\ 



1/2-1 



U2 



ur. 



\B2\ 

\Bi\ 



, (50) 



where transfer matrix T2 is: 



' _ /3J||(k) 
•^ 2 

:__ ;3J||(k) 
2 



4 4 



/3J||(k) 



(51) 

IA2 is matrix built of eigenvectors of T2 and eigenvalues 
A+, A_ of T2and read: 



A± = 



c- 



PJ\\ (k) 



p\Jl 



± 



/3J|| (k) 



Finally, we arrive at the value of |i?„| determinant: 



|i3J = ^fA!+A^U^(Al-A- 



2\B2\-{pJi 



c- 



PJ\\ (k) 



(53) 



- {PJi^r 



with A_)_ and A_ given by Eq. H52|l . Thus, with the 
result in Eq. H53(l , the calculation of the determinant of 
the matrix A^^ in Eq. H48|l and thereby the partition 
function given by Eq. H35|l is completed. 



IV. CRITICAL TEMPERATURE 

Using the partition function in Eq. H35|l with the def- 
initions in Eqs. (|46|1 and H49|) the free energy / is given 
by: 



c 



/ = ^VlnA,' 



iw, -1 



(k) 



2/3iV|| 



IniVii 



(27r/a)' 



■InAi(k) (54) 



while the saddle-point equation for the spherical La- 
grange multiplier becomes: 



/3J||(k)- 



{PJi 



(52) 



g/(C) 
dQ 



= 0. 



(55) 



C=Co 



Explicitly, 



1 

1 = - 

n 



d^k '-^./Djyi+k'-^ 



(27r/a)' C„(k)/D7(k) + /?„(k) 
where we have defined: 



C=Co 
(56) 



Cn (k) = \B, 
I?„(k) 



\B„ 



PJ. 



(3J, 



l-Bi)-5 



l-Bi)-2 



(/3j;)Vn-l| 



(57) 



The in-plane momentum dependence of Ai (k). C„ (k) 
and Dn (k) comes from J|| (k) dependence of |iJn|, see Eq. 
(|53|l . The saddle-point value Co can be determined from 
the condition that in criticality ordcr-paramctcr suscep- 
tibility becomes infinite (see, Eq. ffT^i ): 



C„(k)vAD;:(k) + Z?„(k) 



dCn (k) 

dc 



VAT^k) 



2 dc 



= 0. (58) 



k = 

C = Co 



Closed formulas for the critical temperature can be easily 
obtained only for small values of n. For example, in the 
case of a bilayer system (n = 2) , the respective determi- 
nants read: 



\B2\^ 



c- 



PJ\\ (k) 



|Si| = C- 
\Bo\ = 1, 



2 

PJ\\ (k) 



PJi 



(59) 



and the saddle point equation in Eq. H56|) can be written 
explicitly: 



1 = 



d^k 



C-i/3J||(k) 



{27T/ay 



Co - I (J|| (k) - Jl) 
1 



/32J! 



Co - 1 ( J|| (k) + j; 



with the Lagrange multiplier: 

/3 



/3V^ 



(60) 



Co 



J|l {k^O) + JjL + J^ 



(61) 




Figure 4: (Color online) Critical temperature vs. number of 
Cu-O layers within block for J±/J\\ ~ 0.1 and J^/J\\ = 0.01 
(blue filled circles). Empty circles (green) denote the critical 
temperature dependence Tc (n) = Tc (1)4-0.28 (1 — 1/n), from 
Ref. I23 (factor of 0.28 was chosen to fit the results of the 
present paper) . Diamonds (red) represent results from Monte- 
Carlo (MC) simulations of the classical XY model (Carlson, et 
al. , Ref. [291 . The data of Ref. |23 is scaled by overall factor of 
0.77 to accommodate the finite size effect of MC simulations. 



Finally, the critical temperature of this system (cf. Eq. 
(23) is given by 



Pc 



1 /•+°° ,, ,,, 2(2-e)j|i + j± + / 



(2-OJ||[(2-0^||+./J 
1 



[2(2-e)J|| + (J± + 2Jl)] -Jl 



(62) 



with 



d^k 



r(5 [C — cos (afcx) — cos (afcj,)] (63) 



(27r/a)' 
being density of states of the square lattice given by: 




where & (x) is the unit-step function. For higher values 
of n we can resort to direct numerical evaluation of Eqs. 
(|5H|l and (|55|) in order to compute Tc (n) . 

The dependence of the critical temperature on the 
number of layers n in a block was presented in Fig. 0] 
Phase stiffnesses J|| , J± and Jj^ were chosen to satisfy 



Jii > Jl > J 



(65) 



which is physically reasonable, since the interlayer cou- 
plings J± and J^ are much smaller than in-plane phase 
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j_;/j| 0.1 



Figure 5: (Color online) Critical temperature Tc in) for n = 4 
as a function of c-axis phase stiffnesses Jj_ and J^ . 

stiffness J\\ and inter-layer (intra-block) phase stiffness 

J_L is greater than inter-block coupling J_|_. We found 
that the critical temperature increases monotonically 
with n. This has a simple explanation: for n = 1 (single- 
layer system) with the fixed phase stiffness J\\, the criti- 
cal temperature is determined by J ^ , see Fig. |21 On the 
other hand, in the limit of rt -^- oo, n-layer blocks are of 
infinite size and inter-block coupling J^ is effectively re- 
placed by J_L. Because J±^ > J^, according to the result 
from the previous Section, Tc (n = cxj) > Tc{n = 1). For 
intermediate values on n, one can expect a monotonic 
increase of the Tc (n) with increasing n, as depicted in 
Fig. ^ In this respect, our findings are similar to the 
results from Ref. |23, where the increase of the critical 
temperature of identically doped members of the same 
homologous series due to the interlayer Coulomb interac- 
tion was found to be 



Tc (n) ^ Tc (1) + const x 1 



1 



(66) 



Similar results were obtained by Monte-Carlo 
simulations f2£ where Tc (n) was computed numeri- 
cally by means of the Binder parameter— for systems 
of size up to 24 X 24 X 24. Additionally, we present 
the influence of inter-planar coupling constants J± and 
J_i_ on Tc, see Figs. O and El For a bilayer system 
the role of both interaction constants is identical. For 
systems with n > 3 the influence of the inter-block 
coupling Jj^ is less visible, however J^ = suppresses 
the critical temperature to zero, since for finite n the 
isolated stack of n coupled layers is effectively quasi 
two-dimensional. Therefore, we conclude that the 
change of the number n of layers within block cannot 
alone lead to the appearance of the maximum in the 
critical temperature Tc as a function of n for fixed values 




Figure 6: Critical temperature Tc (n) as a function of c-axis 
phase stiffnesses J± and J^ for various values of n as indicated 
in the inset. 



of the microscopic phase stiffnesses. It is necessary to 
introduce another factor, which acts competitively to the 
effective increasing of inter- layer coupling (for Jj_ > Jj_) 
triggered by increase of ti. As in homologous series of 
superconducting cuprates in-plane lattice parameters 
and distances between Cu-0 planes are constant among 
different members of the family (see, Ref. [l3) . one 
can expect that the change of n will mainly affect J|| 
coupling parameter. Without assuming any specific 
mechanism or introducing a competing (hidden) order, 
one can calculate values of Jii for given J± and Jj^ 
leading to dependence of critical temperature on n, 
which is in agreement with experimental findings. The 
result is presented in Fig. [7|and in Table UTI It is visible 
that for n — 2, 3 depicted dependence changes from 
rising to falling. This coincides with the fact that for 
n = 3 layers lose their equivalence i.e. they split for 
two groups: outer, lying on block surfaces and inner. 
Because, in real materials, within block of Cu-0 planes 
c-axis redistribution of charge is possible, one can expect 
that the nonequivalence of the planes within a block may 
lead to renormalization of the in-plane phase stiffnesses 
of the respective layers. 



V. CHARGE REDISTRIBUTION 

In this section, using the result of the previous para- 
graphs, we present a phenomenological theory for the 
downturn of the Tc (n) with increasing number of layers 
in multilayer materials and compare it against the ob- 
served systematics, providing support to the presented 
theory. The idea is simple: it is well known that a su- 
perconductor is characterized by a global phase stiffness, 
which is proportional to the condensate density. In high- 
Tc cuprates the global phase stiffness is strongly doping 
dependent and this dependence has a universal character 
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n 


T.[K1 


J\\ [meV] 


1 


97 


19.1 


2 


126 


19.1 


3 


135 


18.3 


4 


125 


15.8 


5 


110 


13.1 


6 


97 


11.0 


7 


88 


9.64 



Table II: Superconducting critical temperature Tc (n) of 
HgBa2Can-iCun02n+2+a compounds taken from Ref. [l3 
and the corresponding values of the in-plane stiffnesses Jy 

calculated from Eq. 1561 for J± = 0.261meV and J^ = 
0.000361meV. 
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Figure 7: In-plane phase stiffness Jy , from Table |n] as a 
function of number of layers n. 



as evidenced experimentally by the dome-shaped depen- 
dence of the critical temperature as a function of dop- 
ing. Phrased in the language of the presented model 
with phase variables and microscopic phase stiffnesses, 
it means that the parameters of this model have to be 
changed substantially with doping, since the global phase 
stiffness (i.e. the condensate) builds on microscopic phase 
stiffnesses and phase variables. The dependence of the 
phase stiffnesses on doping 5 obviously cannot be found 
in the presented model, since they are treated as effective 
values to be determined, in principle, from the (yet un- 
available) microscopic theory. However, assuming that 
all of them share the same doping dependence: 



^11 = S{S)J\\ 

j; = g{5)J^, (67) 



one can rewrite the saddle point equation in Eq. (|56|l 
the form: 



Te(,5)=.g(,5)/i(j||,Ji,j;) 



(68) 



where ft, is a function that is independent on doping 6. For 
a given microscopic phase stiffnesses it follows from Eq. 
(|68|1 that the function g (S) determines the overall dop- 
ing dependence of T^. {S). However, the parabohc shape 
of Tc{S) is well-known (see, Ref. Isij) and the supercon- 
ducting "dome" in the T — S phase diagram extends from 
S Ri 0.05 to (5 « 0.25 with the optimal doping Sopt ~ 0.15. 
This allows us to extract the doping dependence of g (S) 
as a quadratic function of i5 in the form: 



,^,\ l-^iS- Soptf for 0.05 <S< 0.025 
■'^ 1 for (5 < 0.05 or (5 > 0.25, 



(69) 



where a = 0.01. According to results of the previous 
section, in the case of strong anisotropy (J|| ^ J±,J^), 
the in-plane phase stiffness Jy has the major influence on 
the critical temperature Tc, which is weakly dependent 
on inter-planar phase stiffnesses (see, Fig. ^. Therefore, 
we limit further discussion to the case, where the inter- 
planar phase stiffnesses as functions of doping J± « Jj_ 
and Jj^ « Jj^ are approximately constant. The value of 
J|| can be determined using the following procedure. Ob- 
viously, single- and double-layer systems cannot exhibit 
charge redistribution: all the planes arc equivalent and 
therefore, the layers can be optimally doped g (Sopt) = 1- 
Using the data for the critical temperature Tc (n) for 
optimally-doped mercuro-cupratesii for n = 1, 2, we can 
determine Jy and the ratio Jj_/J±. Unfortunately, the 
estimates of the values of couplings along c-axis J± and 
Jj^ are difficult to make. However, using the relation be- 
tween zero-temperature c-axis and ab penetration depths 
and bare phase stiffnesses!^ 



1 



A| (T = 0) 

1 
Ai(T = 0) 



^11 



(70) 



we can evaluate J^/Jh to determine the value of the 
weakest coupling J^. In order to do this, we assume 
the largest anisotropy ratio 



A^/A|| - 230 



(71) 



found in cupratesiS In this way all the phase stiffnesses 
can be fixed; we found that Jy = 19.1meV (see. Table IH)! . 
while J± = 0.26meV and J'j_ = 0.00036meV. Ahhough, 
the estimates of c-axis phase stiffnesses in general can 
substantially vary their influence on the critical temper- 
ature is rather weak, as shown in Section II. 

In the multi-layer systems (n > 3) the planes are 
not equivalent any more and charge redistribution within 
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blocks is possiblei^^iS^ Therefore, in-plane phase stiffness 
of every layer i becomes J\\ {5i). Since our model does 
not contain independent values of J|| for different planes 
(which would complicate theoretical considerations con- 
siderably - given the lack of the translational symmetry 
along c-axis), we modify Eq. H67|l and use the average 
over n layers in the n-layer block: 



J\\ 



J\\^J^i^) = '-}Y.9i5^) 



(72) 



However, nuclear magnetic resonance (NMR) measure- 
ments show that charge distribution among inner planes 
is almost constanti^ Therefore, we characterize charge 
distribution by two parameters: we assume that the dop- 
ing of inner and outer planes is given by Sip and 5 op, re- 
spectively, see Fig. |S1 Introducing a ratio of dopings for 
outer and inner planes 



R = 5op/hp > 1 



(73) 



and noting that in every block there are 2 outer and n — 2 
inner layers, we can write Eq. (|72|l as: 

jr (<5,p, i?, 7i) = i [2g (RS.p) + in-2)g (S^p)] . (74) 

By observing that the maximum of the critical 
temperature T^ (n) corresponds to the maximum of 
J«'" (Sip, R, n), we can find the optimal value of Sip ~ S^pt 
using the condition 

dJl^{S,p,R,n) 



dS; 



ip 



0. 



(75) 



Explicitly, using Eq. ifTHl we obtain for S*pf; 

^2 n-2 + 2R 
°P* °P*n-2-)-2i?2' 



(76) 



which allows to find the value of the effective in-plane 
phase stiffness 



J||^jr(n,i?) 



(77) 



as a function of number of layers n and charge imbalance 
ratio R: 



J|f" (n, R) = J|| 



1) 



,5,%(n-2)(i?- , 
an (71 - 2 + 2i?2) 



(78) 



Using Eq. H78(l we can determine the charge imbal- 
ance parameters i? as a function of number of layers n 
for a given set of experimental data for the critical tem- 
perature Tc (n) . For this purpose we use the dependence 
of the superconducting critical temperature on number 
of layers reported by Kuzemskaya. et al.. Ref. [T3- The 
procedure reads as follows: 




/ outer plane 
/ overdoped 



inner planes 
underdoped 



outer plane 
ovsrdoped 



Figure 8: (Color online) Schematic representation of charge 
redistribution leading to overdoping of outer and underdoping 
of inner layers, respectively. 



(i) For n = 1 and n = 2 by definition R = 1 and 
consequently J«-'" {n, R) = Jy. Using the value of Tc (1) 

from experiment, we can determine the in-plane J|| and 
Jj^ = J± (for n = 1) phase stiffnesses using Eq. (|5^ for 
the critical temperature Tc{n = 1) and the supplemen- 
tary equations ifTUIl and (|7T|) . 

(ii) For n = 2 we determine in-block J_l phase stiffness 
using Eq. (|56|l for Tc {n = 2), while keeping J|| and J^ 
unchanged. In this way all the energy setting parameters 
J|| , J^ and J^ are fixed. 

(iii) For n > 2, with the computed parameters Jy, Jj_ 

and Jj^, by making use of Eq. (|56|1 for theoretical Tc (n), 
we calculate the values of R, using Tc {n) as an input 
from experiment. 

The obtained values of the ratio R are presented in Fig. 
|51 while resulting dopings of inner and outer planes are 
shown in Fig. 1101 From these results we can infer that for 
n = 1, 2, 3 the critical temperature is rising due to the in- 
creased inter-layer (intra-block) phase stiffness J± > J^, 
while the charge redistribution for n = 3 is rather small. 
For values of n > 3 the critical temperature is very sen- 
sitive to charge redistribution and we have a decrease 
from Tc = 135K for n = 3 (i? = 1.35) to Tc = 125K 
for n — i {R — 1.82). The calculated values of R are 
in good agreement with experimental results based on 
NMR measurement of HgBa2Ca4Cu50j, compound by 
Kotegawa, Ref. |3J. Interestingly, the obtained parame- 
ters R for n = 3, 4 are much smaller than those found 
in Ref. |3^ [R — 4.5 and i? = 9, respectively) calculated 
using simplified considerations based on competition be- 
tween electrostatic and band energy estimates. However, 
in the real materials charge depletion in the inner layers 
may be high enough to destroy the superconductivity in 
these layers, i.e. doping of one type of planes may be well 
outside superconducting region. This scenario have also 
been confirmed experimentally by NMR studies of five- 
layer compound HgBa2Ca4Cu50yi^ Only outer Cu-0 
planes are superconducting with Tc — 108/'S'. while inner 
planes exhibit a static antiferromagnetic ordering with 
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n 


i? (this work) 


R (Ref. ^ 


R (Ref. ^ 


R (Ref. ^ 


1 


1 


1 


1 


n.a. 


2 


1 


1 


1 


n.a. 


3 


1.35 


1.14 


4.5" 


n.a. 


4 


1.82 


1.49 


9" 


n.a. 


5 


2.23 


1.64 


n.a. 


3.8"^ 


6 


2.56 


n.a. 


n.a. 


n.a. 


7 


2.78 


n.a. 


n.a. 


n.a. 



Table III: Calculated values of the ratio R (this work, second 
column) as a function of number of layers n using the data for 
Tc (n) from Kuzemskaya, et al. Ref. [ITI . compared with the re- 
sults from other works: R based on calculated values of charge 
concentration of outer (0.45"''') and inner layers (0.1", 0.05*') 
using a theory of competition between electrostatic and band 
energy estimates (Di Stasio, et al., Ref. 1351) : nuclear magnetic 
resonance measurement of HgBa2Ca4Cu50y compound'^ with 
the values of charge concentration of outer (0.212) and inner 
layers (0.057); n.a. - non-available. 



0.25 



0.20. 



0.15 



0.10. 



0.05 



♦ imier plaiies 
n oi-iter planes 



4 

n 



Figure 10: Doping of inner Sip and outer Sop planes vs. num- 
ber of layers n that corresponds to the data from Fig. |^ 
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Figure 9: Critical temperature Tc_ vs. number of layers n that 
fits the experimental data from Ref. [13 and the calculated 
values of parameter R describing the ratio of doping of outer 
to inner planes. Dashed line is only guide for the eye. 



Tjv = &QK ■ This suggest possible extension of our model 
to include independent in-plane phase stiffnesses of the 
inner and outer layers. 



VI. SUMMARY AND CONCLUSIONS 

We have considered a model of n-layer high- 
temperature cuprates of homologous series like 
HgBa2Ca„_iCu„02+2?i+(5 to study the dependence 
of the critical temperature T^ (n) on the number n 
of Cu-0 planes in the elementary cell. We focus on 
the "phase only" description of the high-temperature 
superconducting system motivated by the experimental 



evidence that the ordering of the phase degrees of 
freedom is responsible for the emergence of the super- 
conducting state with long-range order. To this end, 
we have proposed a three-dimensional semi microscopic 
XY model with two-component vectors that involve 
phase variables and adjustable parameters represent- 
ing microscopic phase stiffnesses. The model fully 
implements complicated stacked plane structure along 
c-axis to capture the layered composition of homologous 
series. This allows us to go well beyond widely used 
phenomenological Lawrence-Doniach model for layered 
superconductors. Implementing spherical closure rela- 
tion we have solved the phase XY model exactly with the 
help of transfer matrix method and calculated Tc (n) for 
chosen system parameters and arbitrary block size n. In 
this way we were able to elucidate the role of the c-axis 
anisotropy and its influence on the critical temperature. 
Furthermore, by making a physically justified assump- 
tion regarding the doping dependence of the microscopic 
phase stiffnesses we were able to calculate the bell 
shaped curve of Tc (n) as a function of block size n with 
a maximum at n = 3 by effectively accommodating 
inhomogeneous charge distribution among planes. The 
success of our approach suggests to us that at least for 
the superconducting phase transition properties studied, 
one does not need to invoke any competing "hidden" 
order parameter that appear in other approaches^ In 
the present work, the vanishing of the superconductivity 
with underdoping can be understood by the reduction 
of the in-plane phase stiffnesses and can be linked to the 
manifestation of Mott physics (no double occupancy due 
to the large Coulomb on-site repulsion) that leads to the 
loss of long-range phase coherence while moving toward 
half-filled limit {S = 0). For doping S = 0, the fixed 
electron number implies large fluctuations in the conju- 
gate phase variable, which naturally translates into the 
reduction of the microscopic in-plane phase stiffnesses 



14 



and destruction of the superconducting long-range order 
in this Umit. In the opposite region of large 6, the onset 
of a pair-breaking effect (at the pseudogap temperature 
T*) can deplete the microscopic phase stiffnesses, thus 
reducing the critical temperature. Therefore, it would 
be interesting to relate the phase stiffnesses parameters 
used in the present paper with the microscopic material 
characteristics of the electronic system (like hopping 
parameters, antiferromagnetic exchange and Coulomb 
energy) to establish a link between our semi-microscopic 
approach and physics of strongly correlated systems. 
The charge redistribution taking place in multi-layered 
cuprates also points out the profound impact of charging 
effects and unusual role of the internal electric fields in 
cuprate superconductorsi^ 
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Appendix A: MATRIX DIAGONALIZATION 

The transfer matrix method requires raising transfer 
matrix to an arbitrary power. If matrix T of the or- 
der m has m distinct eigenvalues Ai,...,Am it can be 
diagonalized: 



UVW 



(Al) 



where U is a matrix built of eigenvectors ei, ...,£„ of T: 



U 





ei 




im 



(A2) 



while 



V = 



Al 

■■. 

A™ 



(A3) 



is diagonal matrix. In result any integer power p of T 
can be simply found: 

T'P = (UVU-^) ... (UVU^^) ^ UVPU~\ (A4) 



since U ^U = 1. Finally: 



TP =U 



A? 

■•. 

AE, 



W 



(A5) 



This procedure is utilised in Sec. III. The columns of the 
matrix lAi in Eq. 1421 contain eigenvectors of matrix Ti; 



where: 



Ui 



A\B^\ + {[iJ'^ |B„_2|-%/A 4|B„| + f/3Jl) |S„_2| + VA 



8|B„-i| 
1 



8|-B„_i| 
1 



(A6) 



A 



The inverse of the matrix lAi is given by: 



16|S„|+ /3J^ |B„_ 



- 4/3J^|B„_i 



Ui 



4|B„_i| 4|B„| + (/?j;)Vn-2|+\/A 



^-1^ VA 2|B„_i| 

4|B„„i| 4|i?„| + (/3Jl)V«-2|-VA 
\/A 2|S„^i| 

Similarly, the matrix U2 in Eq. |^ is built of eigenvectors of T2 : 



Kc-^)-h (c-^) -(/^^-)^ K^-^) + ^\ (c-^) -(/5^-r 



(A7) 



(A8) 



(A9) 
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